Dynamics of a ferromagnetic domain wall and the Barkhausen effect 
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We derive an equation of motion for the the dynamics of a ferromagnetic domain wall driven 
by an external magnetic field through a disordered medium and we study the associated depirming 
transition. The long-range dipolar interactions set the upper critical dimension to be d c = 3, 
so we suggest that mean-field exponents describe the Barkhausen effect for three-dimensional soft 
ferromagnetic materials. We analyze the scaling of the Barkhausen jumps as a function of the field 
driving rate and the intensity of the demagnetizing field, and find results in quantitative agreement 
with experiments on crystalline and amorphous soft ferromagnetic alloys. 
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The magnetization of a ferromagnet displays discrete 
jumps as the external magnetic field is increased. This 
phenomenon, known as the Barkhausen effect, was first 
observed in 1919 by recording the tickling noise pro- 
duced by the sudden reversal of the Weiss domains |Q. 
The Barkhausen effect has been widely used as a non- 
destructive method to test magnetic materials and the 
statistical properties of the noise have been analyzed in 
detail In particular, it has been observed that the 

distributions of sizes and durations of Barkhausen jumps 
decay as power laws at low applied field rates Q-|6|. In 
addition to its practical and technological applications, 
the Barkhausen effect has recently attracted a growing 
interest as an example of a complex dynamical system 
displaying critical behavior J^-pT[ . 

In soft ferromagnetic materials the magnetization pro- 
cess is composed of two distinct mechanisms |l2| ]. (i) 
When the field is increased from the saturated region, 
domains nucleate in the sample, typically starting from 
the boundaries, (ii) In the central part of the hysteresis 
loop, around the coercive field, the magnetization pro- 
cess is mainly due to domain wall motion. The disorder 
present in the material (due to non magnetic impurities, 
lattice dislocations, residual stresses, etc.) is responsi- 
ble for the jerky motion of the domain walls, giving rise 
to the jumps observed in the magnetization. The mov- 
ing walls are usually parallel to the magnetization (180° 
walls), and span the sample from end to end pl|. The 
statistical properties of the Barkhausen effect are nor- 
mally studied in the central part of the hysteresis loop 
and can therefore be understood by studying domain wall 
motion. 

It has recently been proposed to relate the scaling 
properties of the Barkhausen noise to the critical be- 
havior expected at the dcpinning transition of an elastic 
interface PJTl|]. The numerical values of the scaling ex- 
ponents, however, do not agree with most experimental 
data [||-^) . Interestingly, a quantitative description of the 
phenomenon can be obtained by a simple phenomenolog- 
ical model where the wall is described as a sing le point 
moving in a correlated random pinning field |13|] . 

Here, we present an accurate treatment of magnetic 



interactions in the context of the depinning transition, 
which allows us to explain the experiments and to give a 
microscopic justification for the model of Ref. [ ^3[ . We 
study a single domain wall separating two regions with 
magnetization of constant magnitude M and opposite di- 
rections. We assume that the wall surface does not form 
overhangs, and describe the position of the wall by its 
height h(r,t). The motion of the wall is overdamped 
because of eddy currents so, neglecting thermal fluctua- 
tions, the evolution of h(f,t) is governed by 



dh(r,t) SE({h(r,t)}) 



dt 



Sh{f, t) 



(1) 



where E({h(f,t)}) is the energy of the system p3 , 15[ . 
We will show that incorporating the effects of ferromag- 
netic, magnetocrystalline, magnetostatic and dipolar in- 
teractions, and disorder, the equation of motion is 



dh(r, t) 
dt 



vvV 2 h{r, t) + H-H d 



+ / d 2 r'K(r-f')(h(f f ',t)-h[f,t))+ri(f f ,h). (2) 



Here Uq is the surface tension, H is the external magnetic 
field, the demagnetizing field H c i and the non-local term 
are due to dipolar interactions, and we model the disor- 
der with a random force, with Gaussian distribution and 
short range correlations 



(rj(r, h)r](f, h')) = 5 2 {r- f')A(h - ti), 



(3) 



where A(x) decays very rapidly for large values of the 
argument. 

The dipolar interactions are treated considering effec- 
tive magnetic charges induced by the discontinuities of 
the magnetization across the boundaries of the sample 
and the domain wall ]l6(| . The corresponding magnetic 
surface charge is given by 



a = (Mi - M 2 ) ■ h, 



(4) 



where n is normal to the surface and Mi and Mi are 
the magnetization vectors on each side of the surface. 
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The surface charge induced at the boundary of the sam- 
ple gives rise to a demagnetizing field which opposes the 
external field Jl2| . The simplest approximation is to con- 
sider this field to be constant throughout the sample, and 
to be proportional to the total magnetization pUpf . This 
yields the term Hd = —k f d 2 r h(r,t) in Eq. (g)7 where 
the demagnetizing factor k takes into account the geome- 
try of the domain structure, the shape of the system and 
its size. 

Magnetic surface charges will also appear on the do- 
main wall when its surface is not parallel to the magne- 
tization. In the limit of infinite anisotropy Jl7[ and small 
bending of the surface, we can express the surface charge 
as 



cr(r) = 2Mcos9 ~ 2M 



dh(r,t) 
dx 



(5) 



where 9 is the local angle between the vector 
normal to the surface and the magnetization (see 
Fig. [I]). The interaction energy of these charges — 
(1/2) J d 2 rd 2 r' a{r)a(f')/\r — r'\ gives rise to the non- 
local kernel present in Eq. (0) HMpfl: 



K(r-r') 



2M" 



Z{x-x'f 



(6) 



The interaction (g) is long-range and anisotropic, as can 
also be seen by considering the Fourier transform 



K(p,q) 



2M A 



Vf 2 



(7) 



where p and q are the two components of the Fourier 
vector. Moreover, an estimate of the order of magnitude 
of K(f— f') shows that it dominates over the surface 
tension for all length scales of interest [ft9[ . 

Apart from the non-local kernel, Eq. (g) is the equation 
proposed in Ref. || which, in turn, reduces when k = 
to a driven elastic interface in the presence of quenched 
disorder |2C]-|2l. 

For k = (no demagnetizing field), Eq. (||) displays 
a depinning transition: i.e. there exists a critical field 
H c such that for H < H c the interface is pinned while 
for H > iJ c , it moves with non-zero average velocity. At 
H = H c , the system exhibits scaling properties: the in- 
terface moves by avalanches whose sizes s and durations 
T distributions follow power laws 



POO 



P(T) - T~ 



(8) 



When k > 0, the demagnetizing field provides an addi- 
tional restoring force that keeps the interface at the de- 
pinning transition PJTl[| if the external field is increased 
adiabatically. 

Using the functional renormalization group scheme in- 
troduced in Ref. p3| , we find that — due to the long- 
range kernel K{r — r') in Eq.(^) — the critical behavior of 
Eq. (H) differs from that of elastic interfaces: the upper 



critical dimension becomes d c — 3 |lj,|24|,^5|] , instead of 
d c = 5 pl| , ^2| . Hence we predict that for d = 3, the mo- 
tion of the domain wall will be described by mean-field 
theory (apart from logarithmic corrections), which yields 
]l8| , p2[ | r = 3/2, a = 2, and that the surface is flat (the 
roughness exponent £ is zero). These values differ sig- 
nificantly from the results of elastic interfaces which in 
d = 3 are r ~ 1.3, a ~ 1.5 and C ^ 0.7 ||||. 

Next we make contact between our approach and the 
conventional approach that reduces the domain wall to a 
single point moving in a random pinning field ||,|l3|,^]. 
To this end, we introduce an infinite range version (d — > 
oo) of Eq. (||), which should have the same critical be- 
havior as Eq. (||) but has the advantage of being much 
simpler to analyze. 

To treat the infinite-range model, we discretize the in- 
terface and consider that all N elements are at the same 
distance from each other. Eq. (0) then becomes p8| 



dhjjt) 
dt 



H{t)- X h + J{h-h i {t)) + m {h), 



(9) 



where h = Y,Z=i h i/ N > X = Nk, J = {v Q + 2M 2 ), and 
the external field H{t) increases at a finite constant rate. 
Summing Eq. (^|) over all sites i, we obtain an equation 
for the total magnetization m = Nh 



dm 
~~dt 



N 



ct — x m + 



(10) 



where the time dependence of the field has been made ex- 
plicit. We can approximate ^ r/i by an effective random 
pinning field W(m), depending only on the magnetiza- 
tion. When the interface moves between two configura- 
tions, the change in W is 



W(m') - W{m) =£'Ar?i, 



(11) 



where the sum is restricted to the sites that have moved 
(i.e. their disorder is changed). The total number of such 
sites scales, on average, as \m' — m\, since for d > d c the 
size of an avalanche scales like its area |l8|,^2| . Assuming 
that the Ar]i are uncorrelated and have random signs, we 
find that the effective pinning field is correlated, 



(\W(m') - W{m)\ 2 ) = D\m' - m\, 



(12) 



where D sets the scale of the fluctuations of W . These 
correlations of "Brownian" type have been experimen- 
tally observed in SiFe alloys || . In the model of Alessan- 
dro et al jl3| (ABBM), the domain wall is treated as a 
single point moving in a Brownian correlated field. The 
ABBM model is equivalent to Eqs.([l0|) and (|l^), and pre- 
dicts that the avalanche exponents should depend on the 
field driving rate, with r = 3/2 — c/2 and a = 2 — c 
where c = c/D. In the adiabatic limit c — > 0, we 
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recover the mean-field exponents of the depinning tran- 
sition. Moreover, our results imply that Brownian corre- 
lations in the pinning field do not reflect peculiar long- 
range correlations in the impurities, but represent an ef- 
fective description of the collective motion of the inter- 
face. 

To confirm the results obtained above, we first simu- 
late an automaton version of Eq. (^) in three dimensions, 
applying an adiabatically increasing magnetic field. In 
this limit, we recover, for the distribution of avalanche 
sizes and durations, the results expected for the depin- 
ning transition in mean-field theory |lgf| . 

To overcome the numerical limitations posed by Eq. (|^) 
and obtain an extensive avalanche statistics as a function 
of c and Xi we next simulate the infinite-range model of 
Eq. (||) . We measure the distributions of domain wall ve- 
locities and avalanche sizes and durations (see Figs. || 
and ||, and Jl^]). For c < 1, we find power laws with 
exponents in perfect agreement with the ABBM model. 
The cut-off of the distributions is determined by the de- 
magnetizing fields and diverges for x (T^j. In the 
case of the avalanche size distribution, it scales as x~ 2 
(see Fig. ||). It would be interesting to experimentally 
study the behavior of the cut-off by externally changing 
the demagnetizing field. 

Our results are in agreement with experiments on crys- 
talline (SiFe) HH and amorphous (Co-base and Fe-base) 
ferromagnetic alloys performed at different values of 
c, which yield r = 3/2 — c/2 and a = 2 — c. A direct 
comparison with experiments in which the parameter c 
is not defined, as in the case of a sinusoidal driving field, 
is problematic The result presented in Ref. || (i.e. 

r ~ 1.3) could be due to a finite driving rate (c ~ 0.4). 
In addition, it is important to remark that the present 
theory applies only if domain wall motion is the domi- 
nant magnetization process. This was indeed the case in 
Refs. ^H|^9|], where the noise was recorded only in the 
region of constant permeability around the coercive field. 
A detailed critical discussion of the experimental results 
reported in the literature can be found in Ref. p9fl . 

In particular geometries, typically frames or toroid 
samples, the demagnetizing field is absent (x = 0) ||. 
It is then possible to observe experimentally the depin- 
ning transition of domain walls. It has been reported 
that the average velocity of the domain walls, in differ- 
ent ferromagnetic materials p0[ , increases for H > H c , 
as 

v ~ (H - H c f, (13) 

with (3 = 1, in agreement with the theory (f3 = 1 is ex- 
pected in mean- field theory [p2|). 

We have seen that avalanche distributions can be de- 
scribed by power laws with exponents that do not depend 
on material details. The power spectrum of the noise 
displays instead a more complex structure and does not 
show such a marked universality. At low frequency the 
power spectrum grows with an exponent varying between 



if> = 0.5 for crystalline alloys and ip = 1 for amorphous 
alloys, while at high frequencies it decays with an expo- 
nent varying between^ — —2 for crystals and tp = —1.6 
for amorphous alloys |p|,^|-|6|j29[| . In samples with a single 
domain wall present, the power spectrum was found to 
decay as u~ 2 |J. 

Following the analysis of Tang et al. |3l]] , we obtain at 
low frequency ip = 1 for the power spectrum measured 
on a single site and -0 = when the signal is averaged 
over the whole system. For the averaged spectrum, we 
also find aw -2 decay at large frequencies, due to the 
Brownian properties of the effective pinning field. The 
discrepancies between theory and experiments could be 
due to the presence of many domain walls interacting 
through the demagnetizing field. When a domain wall 
starts to move, the demagnetizing field increases, creat- 
ing a larger pinning force on the other walls. Therefore, 
on short time scales the interactions between the walls is 
irrelevant and should not change the avalanche distribu- 
tions. On larger time scales, this effect may be important 
and could modify the properties of the power spectrum. 
In order to clarify this issue, it would be necessary to 
analyze in detail the dynamics of many coupled domain 
walls. 

The present theory for the Barkhausen effect, based on 
the depinning of a ferromagnetic domain wall, should ap- 
ply to soft ferromagnetic materials, which are frequently 
used in experimental studies of the Barkhausen effect 
[p[-p|,pyr^] . For hard ferromagnets and rare earth ma- 
terials where strong random anisotropics prevent the for- 
mation of domains, disordered spin models could be ap- 
propriate ||. 
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FIG. 1. The domain wall separating two regions of oppo- 
site magnetization. The discontinuities of the normal com- 
ponent of the magnetization across the domain wall generate 
magnetic charges. 
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FIG. 2. The distribution of avalanche sizes in the infinite 
range model as a function of c for N = 32696, \ = 0.0075. 
The lines are the theoretical predictions r = 3/2 — c/2. 
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FIG. 3. The distribution of avalanche sizes in the infi- 
nite range model as a function of \ f° r N = 32696 and 
2 = 0. The data collapse is obtained from the scaling function 
P(s, X ) = s~^f( X 2 s). 
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